Search for Invisible Decays of 77 and rf in J ftp — > 0r/ and ^77' 
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Using a data sample of 58 x 10 6 J/ip decays collected with the BES II detector at the BEPC, 
searches for invisible decays of n and rf in J/ip to (j>n and <j>n' are performed. The <j> signals, which 
are reconstructed in K + K~ final states, are used to tag the n and n' decays. No signals are found 
for the invisible decays of either n or n' , and upper limits at the 90% confidence level are determined 
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to be 1.65 x 1(T 3 for the ratio B(,? R 7 invis ™ c) and 6.69 x 1(T 2 for gt^^^e) _ These are the first 
searches for r\ and r/ decays into invisible final states. 

PACS numbers: 13.25.Gv, 13.25.Jx, 14.40.Aq, 95.30.Cq 



Invisible decays of quarkonium states such as the J/tp 
and the T, etc., offer a window into what may lie beyond 
the Standard Model (SM) HQ. The reason is that apart 
from neutrinos, the Standard Model includes no other 
invisible final particles that these states can decay into. 
It is such a window that we intend to further explore by 
presenting here the first experimental limits on invisible 
decays of the rj and 77', which complement the limit of 
2.7 10 -7 recently established in 3] for the invisible decays 
of the 7T°. 

Theories beyond the SM generally include new physics, 
such as, possibly light dark matter (LDM) particles 0. 
These can have the right relic abundance to constitute 
the nonbaryonic dark matter of the Universe, if they arc 
coupled to the SM through a new light gauge boson U || , 
or exchanges of heavy fermions. It is also possible to con- 
sider a light neutralino with coupling to the SM mediated 
by a light scalar singlet in the next-to- minimal supersym- 
metric standard model Q. 

Recently, observations of a bright 511 keV 7-ray line 
from the galactic bulge have been reported by the SPI 
spectrometer on the INTEGRAL satellite Q. The cor- 
responding galactic positron flux, as well as the smooth 
symmetric morphology of the 511 keV emission, may be 
interpreted as originating from the annihilation of LDM 
particles into e + e~ pairs jj| (also constrained by It 
is in any case very interesting to search for such light in- 
visible particles in collider experiments. CLEO gave an 
upper bound on T(1S) — > 7 + invisible, which is sen- 
sitive to dark matter candidates lighter than about 3 
GeV/c 2 9], and also provides an upper limit on the axial 
coupling of the new U boson to the b quark. It is cru- 
cial, in addition, to search for the invisible decays of light 
quarkonium (qq, q = u,d, or s quark) states which can 
be used to constrain the masses of LDM particles and 
the couplings of the new boson to the light quarks Q. 
We present here measurements of branching fractions of 
77 and 77' decays into invisible final states. 

The data used in this analysis, consisting of 58 x 10 6 
J/tp events, were accumulated with the BES II de- 
tector [lo|. at the BEPC. BES II is a conventional 
solenoidal magnetic detector that is described in detail 
in Ref. 0|- A 12-layer vertex chamber (VC) surround- 
ing the beam pipe provides trigger and coordinate infor- 
mation. A forty-layer main drift chamber (MDC), lo- 
cated radially outside the VC, provides trajectory and 
energy loss (dE/dx) information for charged tracks over 
85% of the total solid angle. The momentum resolution is 
<T p /p = 0.017-^/i + p 2 (p in GeV/c), and the dE/dx reso- 
lution for hadron tracks is ~ 8%. An array of 48 scintilla- 



tion counters surrounding the MDC measures the time- 
of-flight (TOF) of charged tracks with a resolution of 
~ 200 ps for hadrons. Radially outside the TOF system 
is a 12 radiation length, lead-gas barrel shower counter 
(BSC). This measures the energies of electrons and pho- 
tons over ~ 80% of the total solid angle with an energy 
resolution of <r E /E = 22%/ 'yfE [E in GeV). Outside of 
the solenoid coil, which provides a 0.4 Tesla magnetic 
field over the tracking volume, is an iron flux return that 
is instrumented with three double layers of counters that 
identify muons with momentum greater than 0.5 GeV/c. 




FIG. 1: Schematic of J/tp — > <j>rj or cjyq' . The <jf>, which is 
reconstructed in K + K~ final states, can be used to tag the 
invisible decay of the r\ and 77'. 

In order to detect invisible 77 and 77' decays, we use 
J/ij) — > <jnj and (p>v( decays. These two-body decays pro- 
vide a very simple event topology, as shown in Fig. ^ 
in which the <fi signals can be reconstructed easily and 
cleanly decaying into K + K~ . The reconstructed <j> par- 
ticles can be used to tag the 77 and 77' in order to allow a 
search for their invisible decays. Since both <p and 77 (7/) 
have narrow widths, which are negligible compared with 
the detector resolution, the shape of the momentum dis- 
tribution of the <p is approximately Gaussian. The mean 
value of the cf> momentum distribution is 1.320 GeV/c for 
J/tp -> 077 and 1.192 GeV/c for J/tp -> (jnf. The miss- 
ing momentum, P m i ss — \P m iss\i is a powerful discrim- 
inating variable to separate signal events from possible 
backgrounds, in which the missing side is not from 77 (77') 
decay. Here, P m iss — —P$- The 77 and 77' signal regions 
are defined as \P m iss — 1.320| < 3a^ eso for J/tp — > ^77 and 
\P mlss - 1.192| < 3< so for J/il> -> 077', where a? eso (22 
MeV/c) and o]l eso (20 MeV/c) are detector resolutions 
of P m iss for J/tp — » <prj and J/tp — > <pr)', respectively. In 
addition, the 77 and 77' decay regions are easy to define in 
the lab system due to the strong boost of the <p from J / ip 
decay, as shown in Fig. 

In the event selection, the total number of charged 
tracks is required to be two with net charge zero. Each 
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track should have a good helix fit in the MDC, and the 
polar angle 9 must satisfy | cos#| < 0.8. The event must 
originate near the collision point; tracks must satisfy 
\J x 2 + y 2 < 2 cm, \z\ < 20 cm, where x, y, and z are 
the space coordinates of the point of closest approach of 
tracks to the beam axis. Particle identification (PID) is 
performed using combined TOF and dE/dx information, 
and both charged tracks must be identified as kaons. 

We require that events have no other charged tracks 
besides those of the <p — > K + K~ candidate. We count 
the number of BSC clusters, that could be from a K\ 
or a photon, Nbsc, an d require that Nbsc be zero in 
the region outside cones of 30° around the charged kaon 
tracks. These requirements reject most 77 and rj' decays 
into visible final states. They also eliminate most back- 
grounds from multi-body decays of J/ip — ► <f> + anything. 
In order to ensure that 77 and r[ decay particles are in- 
side the fiducial volume of the detector, the recoil di- 
rection against the <p is required to be within the re- 
gion | cos Q reC oii \ < 0.7, where reco u is the polar angle of 
Pmiss- Figure|2a) shows the invariant mass distribution 
of K + K~ candidates, itikk, after the above selection. A 
clear <p peak is seen. Figure HJb) shows the P m iss distri- 
bution for events with 1.005 < tjikk < 1-035 GeV/c 2 . 
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FIG. 2: (a) The itlkk distribution for candidate events. The 
arrows on the plot indicate the signal region of <j> candidates, 
(b) Pmiss distribution for the events with 1.005 < txikk < 
1.035 GeV/c 2 in (a). The means of the missing momenta for 
J/ip — » <pr\ and J/ip — » <j>rj' are located around 1.32 and 1.20 
GeV/c, respectively, as indicated by the two arrows. 



We use Monte Carlo (MC) simulated events to de- 
termine selection efficiencies for the signal channels and 



study possible backgrounds. We obtain efficiencies of 
23.5% and 23.2% for -q and 77' invisible decays, respec- 
tively. More than 20 exclusive decay modes are stud- 
ied with full MC simulations in order to understand 
the backgrounds. The sources of backgrounds are di- 
vided into two classes. Class I: the background is from 
J/ip — > (f>r)(r)'), where cp — > K + K~ and 77(77') decays 
into other modes than the invisible final states. We 
find the expected number of background events from 
this class is negligible for both 77 and 77'. Class II: 
it is mainly from J/ip decays to the final states with- 
out 77 or 77', such as (pK L K L , ^/ (980)(/ (980) -> 
K L K L ), and K*°K L (K*° -> K^). For 77 case, the 
dominated background is from the decay of J/ip — > 
K*°K L (K*° K+n*), while for 77' case, the dominated 
background is from the decays of J/ip — > (pK^K^ and 
(pfo (980) (f (980) -> K L K L ). The expected number of 
background events from class II is 3.0 ± 0.2 and 90 ± 64 
in the signal regions for 77 and ?/, respectively. 

An unbinned extended maximum likelihood (ML) fit 
is used to extract the event yield for J/ip — > <pf](f]') 
\(p -» K+K- and 77(77') -> invisible]. In the ML fit, 
we require that 1.00 < P m iss < 1-45 GeV/c, shown 
in Fig. Gib), where the background shape is well un- 
derstood. We construct probability density functions 
(PDFs) for the P miss distributions for (J^ ig and T^ lg ) 
signals and background (T^kgd) using detailed simula- 
tions of signal and background. The PDFs for signals 
are parameterized by double Gaussian distributions with 
common means, one relative fraction and two distinct 
widths, which are all fixed to the MC simulation. The 
PDF for background is a bifurcated Gaussian plus a first 
order Polynomial (Pi). All parameters related to the 
background shape are floated in the fit to data. The 
PDFs for signals and background are combined in the 
likelihood function C, defined as a function of the free 
and N bkgdl 



parameters N^ ig 



-(N^ g +N\ g + N bkgd ) 

~N~l 

n[JV' T v - (P l ■ ) + 
l ly sig J sigV miss) ' 



N 



^siq-^siniPmiss) + ^bkgd^bkgd (Pmiss )] ' 



(1) 



where N n slg 
K+K 



and N^ ig are the number of J/ip — > </>(—► 
)r)(— > invisible) and J/ip — > (j>(— » K + K~)r]' (—> 
invisible) signal events; N^gd is the number of back- 
ground events. The fixed parameter N is the total 
number of selected events in the fit region, and P miss 
is the value of P m iss for the «th event. The negative 
log- likelihood (— ln£) is then minimized with respect to 



sig ' sig 

events are used in the fit, and the resulting fitted val 



iVlL, and Nbkgd in the data sample. A total of 105 



ues of N? ig , 



sig 1 



and Ni 



bkgd 



-2.8 ± 1.4, 2.2 ± 3.4, 
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and 106 ± 11, where the errors are statistical. Figure [3| 
shows the P m iss distribution and fitted result. No signif- 
icant signal is observed for the invisible decay of either 
X] or 77'. We obtain upper limits by integrating the nor- 
malized likelihood distribution over the positive values of 
the number of signal events. The upper limits at the 90% 
confidence level are 3.56 events for 77 and 5.72 events for 
rj , respectively. 
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FIG. 3: The Pmiss distribution for candidate events. The 
data (black crosses) are compared to the total fit results. The 
dotted curve is the projection of r\ signal component, and the 
dashed curve is the the projection of rj signal component, and 
the solid curve is the total likelihood fit result. 

The branching fraction of 77(77') — ► 77 is also deter- 
mined in J/ifi — » <prj{rj') decays, in order to obtain the 
ratio of 6(77(77') — > invisible) to 73(77(77') — > 77). The ad- 



vantage of measuring 



6(77(77') — > invisible) 



is that the un- 



6(77(77') -> 77) 
certainties due to the total number of J/ifi events, track- 
ing efficiency, PID, the number of the charged tracks, the 
cut on rriKK-, and residual noise in the BSC cancel. 

The selection criteria for the charged tracks are the 
same as those for J /if) — > ^77(77'), 77(77') — + invisible de- 
cays. However, at least two good photons are required. 
A candidate photon must have hits in the BSC. The num- 
ber of layers hit must be greater than one, and the de- 
posited energy in the BSC more than 50 MeV. The angle 
between the photon emission direction and the shower 
development direction of the neutral track in BSC is re- 
quired to be less than 25°. The opening angles between 
the candidate photons and the charged tracks must be 
greater than 30°. 

The events are kinematically fitted using energy and 
momentum conservation constraints (4-C) under the 
J/ifi — > KKjj hypothesis in order to obtain better mass 
resolution and suppress backgrounds further. We require 
the kinematic fit X 2 K + K - 11 less that 50 (15) for the 77 (77') 
case. If there are more than two photons, the fit is re- 
peated using all permutations, and the combination with 
the best fit to K + K~jj is retained. The numbers of 
J/ifi — > 4>r)(r)')[(j) — ^ K + K~ and 77(77') — > 77] events are 
obtained from fits to the 77 invariant mass distributions. 



The fitted results for 77(77') — > 77 are shown in Fig.QJ 

Contributions to the systematic error on the ratios are 
summarized in Table [fl Systematic errors in the ML fit 
originate from the limited number of events in the data 
sample and from uncertainties in the PDF parameteriza- 
tions. The uncertainty due to the background shape has 
been estimated by varying the PDF shape of the back- 
ground in the ML fit. 

The uncertainty, due to the requirement of no neutral 
clusters in the BSC allowed outside the 30° cones around 
the charged tracks, is obtained using the control sample 
of fully reconstructed J/ifi — > (fir], rj — > 77 events. The 
ratios of events with the requirement on the number of 
extra photons to events without the requirement are ob- 
tained for both data and MC simulation. The difference, 
5%, is considered as the systematic error for both the 77 
and 77' cases. This study determines the difference in the 
noise in the BSC for MC simulation and data. Compared 
with 77 — * invisible decay, we expect that more noise is 
introduced by the photons in 77 — > 77 decay. So it is 
a conservative estimation of the systematic error due to 
the requirement of no clusters in the BSC for the invisible 
decays of 77 and 77' . 

The uncertainty in the determination of the number of 
observed J/tfi — ► (fii](r]'), <fi — ► K + K~ , 77(77') — * 77 events 
is also estimated. Different background shapes are tried 
in the fit to the 77 invariant mass, and the variation of 
the fitted yields is regarded as a systematic error, which 
is 2.0% (1.0%) for the 77 (77') case. The relative systematic 
bytb 

is about 4.0% [ll|. 

constraint is estimated to be 1.0% (5.2%) |12( for the 77 
(77') case. The uncertainty from the trigger efficiency is 
also considered. The total systematic error, a^ ys (cr^ s ), 
on the ratio is 7.7% (11.1%) for 77 (77'), as summarized in 
Table D 

The upper limit on the ratio of the 6(77 — > invisible) to 
6(77 — » 77) is calculated with 



error caused by the uncertainty of the photon efficiency 
The uncertainty due to the Xk+k-j-j 



B(r) — ► invisible) 
6(77 -> 77) 



< 



n UL/ e V 

m 1 



A77 (l-ft?) 



(2) 





0.20 0.30 0.40 0.50 0.60 0.70 
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FIG. 4: (a) The fit of the 77 invariant mass distribution of 
J/ifi — * (fir], rj — > 77. The dashed line shows the background, 
and the solid line is the total fit result, (b) The same plot but 
for J/ifi — » (fir]' , rj' — » 77. 
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TABLE I: Summary of relative systematic errors. The first 
three lines are for J/tp — » tj>r)(r]'), 77(77') — > invisible. The next 
three are for J/ip — > (j>ri(rj'), r)(r)') — * 77. 



Source of Uncertainties 


sys. error (%) 

v n 


PDF shapes in the ML fit 


3.4 


7.3 


MC statistics 


f.O 


1.0 


Requirement on Nbsc 


5.0 


5.0 


Photon efficiency 


4.0 


4.0 


4-C fit for 77(77') — > 77 


1.0 


5.2 


Background shape for 77(77') — > 77 


2.0 


1.0 


Total 


7.7 


11.1 



where is the 90% upper limit of the number of 

observed events for J/i/j — > (f>r], <f> — > K + K~, 7/ — > 
invisible decay, is the MC determined efficiency for 
the signal channel, n 77 is the number of events for the 
J /if) — > 0?7, — * K + K~, 77 — ► 77 decay, e 77 is the 
MC determined efficiency for the decay mode, and er,, is 

ct^ s ) 2 + (cr,f at ) 2 = 8.1%, where a* y s and of Qt arc the 
total relative systematical error for the 77 case from Ta- 
bleHJand the relative statistical error of n!J 7 , respectively. 

For 77', <v is \/«H 2 + K* Qt ) 2 = 21.6%. The rela- 



tive statistical error of the fitted yield for J/tp — > 077(77'), 
77(77') — > 77, is 2.8% (18.5%) according to the results from 
the fit to the invariant mass of 77 in Fig. We also ob- 
tain the upper limit on the ratio of the B(r)' — > invisible) 
to B(j]' — > 77) by replacing 77 with 77' in Eq. PJ. Since 
only the statistical error is considered when we obtain 
the 90% upper limit of the number of events, to be con- 
servative, njj L and n{j L are shifted up by one sigma of 
the additional uncertainties (a v or ay). 



TABLE II: The numbers used in the calculations of the ra- 
tios in Eq. @, where (rip L ) is the upper limit of the 
signal events at the 90% confidence level, e v (£„/) is the se- 
lection efficiency, n 77 (n!L) is the number of the events of 

J/rl> -» Mv'), <t> -» K+K-, 77(7?') - 77, 4y (e£) is its se- 
lection efficiency, a^ tat (a^ t , at ) is the relative statistical error 

of n 77 (n!L) and cr v (<?y) is the total relative error. 



quantity 


value 
77 77' 


«'j 

fe 17 ' ) 

c 77 V c 77y 
stat / stat\ 

a v (a n , ) 
ff, (cry) 


3.56 5.72 
23.5% 23.2% 
1760.2 ± 49.3 71.6 ± 13.2 
17.6% 15.2% 
2.8% 18.5% 
8.1% 21.6% 



Using the numbers in Table [HJ the upper limit on 
the ratio of 73(77(77') — > invisible) and 73(7/(77') — > 77) 
is obtained at the 90% confidence level of 1.65 x 10~ 3 
(6.69 x 10~ 2 ). 



In summary, we search for the invisible decay modes 
of 77 and 77' for the first time in J/t/j — > <jrr]{rf) using 
the 58 million J/tp events at BES II. We find no sig- 
nal yields for the invisible decays of 77 and 77', and ob- 
£> (77(77') — > invisible) 



tain limits on the ratio, 



The upper 



6(77(77') -> 77) 
limits at the 90% confidence level are 1.65 x 10 -3 and 
73(77 — * invisible) B(rf — > invisible) 



6.69 x 10~ 2 for 



73(77 ^ 77) _ B(t)' 77) _ 

respectively. The advantage of measuring the ratios in- 
stead of the branching fractions of the invisible decays is 
that many uncertainties cancel. 
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